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Abstract 

We study the light quark mass dependence of the D and D s meson decay constants, fjj and fjj s , using a 
covariant formulation of chiral perturbation theory (xPT) at next-to-next-to-leading order (NNLO). Using 
the HPQCD lattice results for the D(D S ) decay constants as a benchmark we show that covariant xPT can 
describe the HPQCD results better than heavy meson ^PT (HMxPT) at both NLO and NNLO. Within the 
same framework, taking into account sub-leading (1/rag, with tjiq the heavy quark mass) corrections to 
the values of the low-energy constants and employing the lattice QCD results for gBB'-w, we estimate the 
ratio of f B JfB to be 1.22^^, which agrees well with the HPQCD result 1.226(26). 

PACS numbers: 12.39.Fe, 13.20.Fc, 14.40.Lb, 12.38.Gc 
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I. INTRODUCTION 



The decay constants of charged pseudoscalar mesons 7r ± , K ± , D ± , Df and B ± play an impor- 
tant role in our understanding of strong interaction physics, e.g., in measurements of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements and in the search for signals of physics beyond the 
standard model (SM). At lowest order, the decay width of a charged pseudoscalar P ± with valence 
quark content qiq 2 decaying into a charged lepton pair (i^vt) via a virtual W ± meson is given by 



T(P± _> e„ e ) = ^-f P m\M P \\ - j^j \V qiq2 \\ (1) 

where mi is the ^ mass, \ V qiq2 | is the CKM matrix element between the constituent quarks qiq 2 
in P ± , and Gp is the Fermi constant. The parameter f P is the decay constant, related to the wave 
function overlap of the q\q 2 pair. Measurements of purely leptonic decay branching fractions and 
lifetimes allow an experimental determination of the product \ V qiq2 fp\. A good knowledge of the 
value of either \ V qiq2 | or f P can then be used to determine the value of the other. 

These decay constants can be accessed both experimentally and through lattice Quantum Chro- 
modynamics (1QCD) simulations. While for f n , f K , f D , experimental measurements agree well 
with lattice QCD calculations, a discrepancy is seen for the value of fp, s '■ The 2008 PDG average 
for f Ds is 273 ± 10 MeV yj], about 3a larger than the most precise N f = 2 + 1 1QCD result from 
the HPQCD/UKQCD collaboration [2], 241 ± 3 MeV. On the other hand, experiments and 1QCD 
calculations agree very well with each other on the value of fp>, /o(expt) = 205.8 ± 8.9 MeV and 
/d(1QCD) = 207 ± 4 MeV. The discrepancy concerning fp, s is quite puzzling because whatever 
systematic errors have affected the 1QCD calculation of fp>, they should also be expected for the 
calculation of fp> s . In this context, constraints imposed by this discrepancy on new physics were 
seriously discussed (see, e.g., Ref. [3D). 

However, the situation has changed recently. With the new (updated) data from CLEO [0- 



60 and Babar QVU , together with the Belle measurement (80, the latest PDG average is fp> s = 
257.5 ±6.1 MeV [9jlJ. The discrepancy is reduced to 2.4a. Lately the HPQCD collaboration has 
also updated its study of the D s decay constant [11]. By including additional results at smaller 
lattice spacing along with improved determinations of the lattice spacing and improved tuning o ? 
the charm and strange quark masses, a new value for the D s decay constant has been reported^: 
1 The October 2010 average from the Heavy Flavor Averaging Group (HFAG) is similar: fn s = 257.3 ± 5.3 MeV 



HO] 



A slightly different but less precise value of fp, s = 250.2 ± 3.6 MeV was obtained in Ref. 11211 as a byproduct from 



the study of the D — > K,lv semileptonic decay scalar form factor by the same collaboration. 
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f Ds = 248.0 ± 2.5 MeV. With the updated results from both the experimental side and the HPQCD 
collaboration, the window for possible new physics in this quantity is significantly reduced 111 in . 

An important part of the uncertainties in heavy quark 1QCD simulations comes from chiral 
extrapolations that are needed in order to extrapolate 1QCD simulations, performed with larger- 
than-physical light quark masses, down to the physical point. Recent 1QCD studies of the D 
(D s ) decay constants, both for Nf = 2 + 1 |2j,[l3j] and Nj = 2 [14], have adopted the one-loop 
heavy-meson chiral perturbation theory (HMxPT) (including its partially-quenched and staggered 
counterparts) to perform chiral extrapolations. In particular, the HPQCD collaboration has used 
the standard continuum chiral expansions through first order but augmented by second- and third- 
order polynomial terms in x q = 5 Wq/8(7r/ 7r ) 2 where B = m^/(m u + m d ) to leading order in 
XPT, arguing that the polynomial terms are required by the precision of the data. It is clear that 
the NLO HMyPT alone fails to describe its data. 

HM^PT Ill5l-ll7ll has been widely employed not only in extrapolating 1QCD simulations but 
also in phenomenology studies and has been remarkably successful over the decades (see Ref. Ill8ll 
for a partial review of early applications). In Ref. tlSD, we have argued that a co variant for- 
mulation of %PT may be a better choice for studying heavy-meson phenomenology and 1QCD 
simulations. This was based on the observation that the counterpart in the SU(3) baryon sector, 
heavy baryon %PT, converges very slowly and often fails to describe both phenomenology and 
lattice data (particularly the latter), e.g., in the description of the lattice data for the masses of the 
lowest-lying baryons Il20l l2lll . On the other hand, covariant baryon vPT was shown to provide 
a much improved description of the same data [|22]. Indeed, in Ref. [19] we have shown that for 
the scattering lengths of light pseudoscalar mesons interacting with D mesons, recoil corrections 
are non-negligible. Given the important role played by fjj (Jd s ) m our understanding of strong- 
interaction physics and the importance of chiral extrapolations in 1QCD simulations, it is timely to 
examine how covariant %PT works in conjunction with the HPQCD fc (JdJ data. 

In this letter we study the light quark mass dependence of the HPQCD f D and f Ds results [2^ 
using a covariant formulation of xPT. It is not our purpose to reanalyze the raw 1QCD data because 
the HPQCD collaboration has performed a comprehensive study. Repeating such a process using 
a different formulation of xPT will not likely yield any significantly different results. Instead, 
we will focus on their final results in the continuum limit as a function of m g /m s , with m q the 
3 Although the HPQCD collaboration has updated its study of the fjj a decay constant, it has not done the same for 
the fjj decay constant, and therefore its fu s / Jo ratio remains the same but with a slightly larger uncertainty. For 
our purposes, it is enough to study the HPQCD 2007 data 
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average of up and down quark masses and m s the strange quark mass. These results can be treated 
as quasi-original lattice data because, for chiral extrapolations, the HPQCD collaboration has used 
the NLO HMxPT result plus two polynomials of higher chiral order. Therefore, any inadequacy 
of the NLO HMxPT should have been remedied by fine-tuning the two polynomials. Accordingly 
the extrapolations should be reliable, apart from the fact that the connection with an order-by-order 
%PT analysis is lost. Our present work tries to close this gap. Using the HPQCD continuum limits 
as a benchmark instead of the raw data not only greatly simplifies our analysis but also highlights 
the most important point we wish to make, namely that the covariant formulation of %PT is more 
suitable for chiral extrapolations of 1QCD data than the HMxPT, at least in the present case. 

This paper is organized as follows. In Section II, we introduce the relevant effective chiral 
Lagrangians and calculate the Feynman diagrams contributing to the D(D S ) decay constants up to 
NNLO. In Section III, we show the numerical results and compare them with those of the HMxPT. 
We also estimate the ratio of f B J f B using the values of the low-energy constants (LECs) fixed in 
the present study and employing the lattice results for gsB'-K- A short summary follows in Section 
IV. 



II. THEORETICAL FRAMEWORK 

The decay constants of heavy-light pseudoscalar and vector mesons with quark content qQ, 
with q one of the u, d, and s quarks and Q either the c or b quark, are defined by 

(o\q^QmpM) = -*fpy, (2) 

(0\q^Q(0)\P;(p,e)} = f P *e^ (3) 
where P q denotes a pseudoscalar meson and P* a vector meson. In this convention, fp has mass 



dimension one and f P * has mass dimension two [23]. From now on, we concentrate on the charm 



sector, D, D s , D* , and D*. The formalism can easily be extended to the bottom sector. 

The coupling of the D (D s ) mesons to the vacuum or to Nambu-Goldstone bosons through the 
left-handed current is described by the following leading chiral order Lagrangian: 

41ce = «((c'P;-^)^), (4) 

where a is a normalization constant with mass dimension two, P = (D°, D + , Df), P* = 
(D*°,D* + ,D* S ), mp is the characteristic mass of the P triplet introduced to conserve heavy 



quark spin symmetry in the mg — > oo limit, i.e., rhp, at NLO and mp, at NNLO (see Table 1), and 
u 2 = U = exp[^] with $ the pseudoscalar octet matrix and F their decay constant in the chiral 
limit. We have introduced a dimensionless coefficient d to distinguish the vector and pseudoscalar 
fields, which is 1 if heavy quark spin symmetry is exact. We need to stress that in our covariant 
formulation of ^PT we do not keep track of explicit 1 / vtlq corrections that break heavy quark spin 
and flavor symmetry, instead we focus on SU(3) breaking. This implies that different couplings 
have to be used for D(D*) and B(B*) mesons. In the present work we only need to make such a 
differentiation in calculating diagram Fig. (Id). In Eq. © we have therefore explicitly pointed out 
that d may be different from 1. In all the other places, we will simply set d equal to 1. 

The leading-order (LO) SU(3) breaking of the D meson decay constants is described by the 
following next-to-leading order (NLO) chiral Lagrangians 



£ (3) 



16tt 2 F 2 



b D {(p; 



BP BP 



mp 



ITLp 



(5) 



where b D and b A are two LECs and M = diag(m 2 r , m\, 2m\ 



mt 



Here and in the following 



(• • • ) always denotes the trace in the corresponding flavor space. 

To study the NLO SU(3) breaking, one has to take into account the DD* (D S D*) and DD S 
(D*D*) mass splittings. Experimentally the DD* and D S D* splittings are similar: 



*DD* 



141.4 MeV and A D D , = 143.8 MeV. 



(6) 



Therefore in our calculation we will take an average of these two splittings, i.e., A = (A DD * + 
A DsD *)/2 = 142.6 MeV. It should be noted that the DD* mass splitting is of sub-leading order 
in the 1/rriQ expansion of heavy quark effective theory. The numbers above show that SU(3) 
breaking of this quantity is less than 2%. The mass splitting in principle can also depend on the 
light quark masses but we expect that the dependence of this "hyperfine" splitting should be much 
weaker than that of the D masJ3, mp,, which we discuss below. 

At NLO, the following Lagrangian is responsible for generating SU(3) breaking between the 
D and D s masses Jl9 1: 

£( 2 > = -2c (PP t )(x+) + 2 Cl (P X+ P j ), (7) 



which yields 



m D = m l + 4c (m 2 + 2m 2 K ) - \c\m\ , 



(8) 



This seems to be supported by quenched 1QCD calculations, see, e.g., Refs. 124 



251 
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™ 2 d s = m l + 4c o {ml + 2m 2 K )+4c 1 {ml-2m 2 K ). (9) 

One may implement this mass splitting in two different ways by either using the HPQCD contin- 
uum limits on the D and D s masses Q2Q to fix the three LECs: m , cq, and C\, or taking into account 
only the DD S mass splitting 

- 8ci(m 2 K - ml) = (m 2 Ds - m 2 D + mf,. - m 2 D *)/2 = A s (m D + m Dg + m D * + m D *)/2, (10) 

where we have introduced A s = m Ds — m D « m D * — m D * w (m Ds — m D + m D * — m£>*)/2. 
In the second approach, using the experimental data for m D , m Da , m D *, and m D *, one obtains 
ci = —0.225. We found that the HPQCD continuum limits on the D and D* masses can be 
described very well using Eqs. (18191) . We also found that using Eqs. (|8I9I ) or Eq. (flOl) gives very 
similar results in our analysis of the D (D s ) decay constants. The results shown below are obtained 
using Eq. (fTOl) to implement the SU(3) breaking and light quark mass evolution of the D (D s ) 
masses. 

In order to calculate loop diagrams contributing to the decay constants one needs to know the 
coupling, (?dd*0, with <p denoting a Nambu-Goldstone boson. This is provided at the leading chiral 



order by the following Lagrangian dl9J] 



= igm D (P*u^P - Pu»P* ] ) (11) 

where = i{v)d^u — ud^v)) and m D = 1972.1 MeV, the average of D, D s , D*, and D* 
masses. The coupling g can be determined from the D* + —s- D°n + decay width, which yields 
9dd*-k = 0.60 ± 0.07 II 1911 . At the chiral order we are working, one can take gr,D*<t> = 9dd*-k- 
If heavy quark flavor symmetry is exact, we expect gsB*-K = 9dd*-k- Otherwise deviations are 
expected. We will come back to this later. 

pi 

TABLE I. Numerical values of (isospin-averaged) masses [9] and decay constants (in units of MeV) used 
in the present study. The eta meson mass is calculated using the Gell-Mann-Okubo mass relation: m 2 = 

(\m\ - ml)/ 3. 

mz> mo A s A tub A s (B) A(_B) m n uik m v f n Fo 
1972.1 1867.2 102.5 142.6 5279.3 88.7 47.5 138.0 495.6 566.7 92.4 1.15/^ 
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(a) 



(b) 



(c) 



(A) 



(e) 



FIG. 1. Feynman diagrams contributing to the calculation of fu and fjj g up to next-to-leading order (NLO): 
(a) and (b) are LO and NLO tree level diagrams, loop diagrams (c), (d), and (e) contribute at NLO. 

Up to NLO, the D(D S ) decay constants receive contributions from the Feynman diagrams 
shown in Fig, [p Studies of these decay constants within the framework of HM^PT have a long 
history Il26l - l28ll . Here we are going to present the first covariant %PT calculation. Insertion of the 
mass splittings between D, D s , D*, and D* in the loop diagrams shown in Fig. Q] generates the 
NNLO contributions which are implemented in the present case by making the following replace- 
ments in the NLO results: 



mjj — y mz) + A s , rriD* — > + A, and m^* — > m^ + A + A t 



(12) 



with the values of these quantities given in Table I. It should be noted that there is no new counter- 
term appearing at 0(p A ). 

Computation of the tree-level diagrams Figs, (la, lb) is trivial. Fig. (la) gives a = a/m p with 
mass dimension one for both D and D s . Fig. (lb) yields 

1 



51 — a 

5 2 = a 



16tt 2 F 2 



16vr 2 F 2 



(b A (2m 2 K + ml) + b D ml) 



(b A (2m 2 K + ml) + b D (2m 2 K - ml)) 



(13) 
(14) 



where 5\ is for D and S 2 for D s . 

Diagram Fig. (lc) is the wave function renormalization, from which one can calculate the wave 
function renormalization constants, which can be written as 



dp, 



(15) 



where pi denotes the four-momentum of D (D s ), mi the mass of D (D s ), mj the mass of D* (D*), 
and m k the mass of it, r/, and K. The coefficients are given in Table II. The function (j) w is 
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defined a 



(-2m 2 M (p 2 + m 2 v ) + (m 2 v - p 2 ) 2 + m 4 M ) B (p 



■ i ,m M ,m v 



±F 2 m v . 

+A (my) (-p 2 + m 2 M - m 2 y) + A (m 2 M ) (-p 2 + 3m 2 M + m 2 v ) 



(16) 



where the functions A and B are defined in the Appendix. 

Diagram Fig. (Id) provides current renormalization, which has the following form 

d = ad &,j,k(t>c( m h m 2 j} m 2 k ), 



(17) 



where are given in Table II with i running over D and D s , j over D* and D*, and k over tt, 
r], K. The function <\> c is defined as 

(gm D )m P 



0c(m i ,m y ,m M ) 



(m^- - my) (m 2 - m 2 M + m 2 v ) B (0, m 2 M , m 2 v ) - 2m 2 A (m 2 M ) 



+ (-2m 2 (m 2 M + m 2 v ) + m\ + (m 2 M - m 2 v ) 2 ) B ( 



m?,ml, /// 



, ■ III u . Illy) 



(18) 



It should be noted that Cj vanishes in NLO HM%PT but plays an important role in covariant %PT 
Diagram Fig. (le) also provides current correction 



Ti = a C;, ; .M"'j)//-; 

j=7T,rj,K 



• 



(19) 



with given in Table III. 

TABLE II. Coefficients, £ij t k, appearing in Eqs. (1151171 ). 
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5 To be consistent, the product grriD is only appropriate for NLO. At NNLO, it has to be replaced by g'rriD with 
g' = grriD /m d ~ 0.63 before performing expansions in terms of 1 /m d either to obtain the HM^PT results or to 
remove the power-counting-breaking pieces. The same applies to the calculation of Ci [see Eq. (TTTb l. 
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TABLE III. Coefficients, Qj, appearing in Eq. (fT9l >. 





7T T) K 


D 

D s 


3 1 1 
8 24 4 

-i -i 
u 6 2 



The total results are then 

ft = a(l + 2i/2) + + T< + C f . (20) 

Because of the large Z) meson masses, C» and contain so-called power-counting-breaking 
(PCB) terms. As explained in detail in Ref. [19] one can simply expand these functions in terms 
of 1/rfiD at NLO or l/mjj at NNLO and then remove the PCB pieces. This procedure is in fact 



the same as the extended-on-mass- she 



baryon chiral perturbation theory 11291 



1 (EOMS) scheme. This scheme was first developed for 



3011 and has been shown to be superior to heavy baryon 



XPT in a number of cases, see, e.g., Refs. J22 H31I I32I1. With the full results for Cj and Z { given 
above the expansion can easily be performed and then one obtains C{ and Z if which have a proper 
power-counting as prescribed in Ref. [19]. At the end one finds 

fi = a(l + Z i /2) + 6 i + T i + C i , (21) 

the expression that is used in the actual calculations. By expanding Z { and Cj in terms of l/rhp at 
NLO or l/m D at NNLO and keeping the lowest order in l/m D (1/mo) one can easily obtain the 
corresponding HMxPT results. 



HI. RESULTS AND DISCUSSION 



Before presenting the numerical results, we should make it clear that in our present formulation 
of x?T we have focused on SU(3) breaking in the context of the chiral expansions but we have not 
utilized explicitly heavy quark symmetry that relates the couplings of the D mesons with those of 
the D*, B, and B* mesons. 

In the present case, we encounter three LECs: a, bo, and 6a- At this point, light quark mass 
dependent 1QCD results are extremely useful. By a least-squares fit to the HPQCD results , one 
can fix those three LECs appearing in our calculation. 
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First we treat the D, D s , D*, and D* mesons as degenerate, i.e., we work up to NLO. The 
corresponding results are shown in Fig. 2, where the HM^PT results are obtained by expanding 
our covariant results in terms of l/rhp and keeping only the lowest-order terms. It is clear that the 
co variant results (with x 2 =41) are in much better agreement with the HPQCD continuum limits 
than the HM%PT results (with x 2 — 201) \ This is not surprising because as we mentioned 
earlier the HPQCD collaboration has added second and third order polynomial terms in x q to 
perform their extrapolation. Furthermore one can notice that at larger light quark masses the 
difference between the covariant %PT and the HM%PT results becomes larger. This highlights the 
importance of using a covariant formulation of ^PT in order to make chiral extrapolations if lattice 
simulations are performed with relatively large light quark masses. Similar conclusions have been 
reached in studying the light quark mass dependence of the lowest-lying octet and decuplet baryon 
masses I122ll . 

Taking into account the mass splittings between D, D s , D*, and D* as prescribed by Eq. (fT2l) 
one obtains the NNLO %PT results. Fitting them to the HPQCD extrapolations, one finds the 
results shown in Fig. 3. Compared to Fig. 2, it is clear that the agreement between the covariant 
XVT results with the HPQCD extrapolations becomes even better. Furthermore the covariant %PT 
results (with x 2 — 16) is still visibly better than the HMxPT results (with x 2 = 59), but now the 
difference between the covariant and the HM %PT results becomes smaller. The three LECs in the 
NNLO covariant xPT have the following values: a = 208 MeV, b D = 0.318, b A = 0.166. 

If we had fitted the HPQCD extrapolations by neglecting the loop contributions, we would have 
obtained a even better agreement (x 2 = 9). In Ref. (12211 we also found that the lattice baryon mass 
data could be fitted better with the LO (linear in m q ) chiral extrapolation. But there we found 
that the NLO chiral results in fact describe the experimental data better than the LO (linear) chiral 
extrapolation. This just shows that the lattice baryon mass data behave more linearly as a function 
of light quark masses at large light quark masses and chiral logarithms play a more relevant role 
at smaller light quark masses, as one naively expects. 

Another way of understanding the importance of chiral logarithms is to perform separate fits for 
lattice simulations obtained at different light quark masses. One expects that at smaller light quark 
masses (e.g., < 300 MeV) covariant %PT and HM ^PT results should perform more or less 
similarly. On the other hand, if the light quark masses are larger, covariant %PT should be a better 
6 It should be noted that the absolute value of x 2 as defined here does not have a clear-cut physical meaning. It only 
reflects to what extent the chiral results agree with the HPQCD extrapolations. 
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FIG. 2. (Color online) Light quark mass dependence of fjj (solid lines) and fjj a (dashed lines). The black 
lines show the results of the NLO HMxPT and the blue lines the results of the covariant NLO ChPT. The 
red lines are the continuum extrapolations of the HPQCD collaboration [2]. The ratio r = m q /m s is related 
to the pseudoscalar meson masses at leading chiral order through = 2B§m s r and m 2 K = Bom s (r + 1) 
with Bo = m^/(m u + m^) and m q = (m u + m<i)/2, where m u , m^, and m s are the physical up, down, 
and strange quark mass. 

choice. In Fig. 4, we show the fitted results obtained from fitting the HPQCD extrapolations in two 
different regions of light quark masses, m q /m s < 0.2 (left panel) and m q /m s > 0.2 (right panel). 
It is clearly seen that fitting lattice data with large light quark masses using the HMxPT results 
may give unreliable extrapolations. Here we have used the NNLO HMxPT and covariant ^PT 
results for comparison. The difference will become even larger if the NLO x?T results are used. 
We should also mention that even for m q /m s < 0.2 (m w < 307 MeV) the HPQCD extrapolations 
are better described by covariant %PT than by HMxPT judging from the x 2 analysis (although the 
difference is so small that it can hardly be appreciated by just looking at the left panel of Fig.©. 

We have checked that our covariant results are stable with respect to variations of certain input 
parameters within reasonable ranges, e.g., m p < /j, < 2 GeV and 0.53 < g < 0.67, where /x is 
the renormalization scale and g the DD*% coupling defined in Eq. (ITT]) . With our standard choice: 
g = 0.6 and /i = 1 GeV, we have also noticed that for the NNLO covariant %PT to produce a 
smaller \ 2 man me linear chiral extrapolation, d has to be larger than 1.23. If we use the quenched 
1QCD result, d = 1.35 ± 0.06 [33], the fit is even bettej^. On the other hand, our results remain 



7 Using the results from a more recent calculation by the UKQCD collaboration Il34ll . one obtains c' f=s 1.18 ± 0.13 
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FIG. 3. (Color online) Same as Fig. 12 but the chiral expansions are calculated up to NNLO. 

qualitatively the same with either d = 1 or d = 1.35. Therefore we have presented the results 
obtained with d — 1 . 

Chiral perturbation theory not only helps extrapolating 1QCD simulations to the physical light 
quark masses. It also benefits from this process because once the values of the relevant LECs are 
fixed by fitting the 1QCD data, ;\PT predicts observables involving the same set of LECs. In the 
present case, assuming that the 1/mg corrections to the values of the three LECs bp, 6^, and g are 
small, we can calculate the ratio of /b s //b by making the following replacements in our NNLO 
covariant %PT results: 



m D ->m B , A^A(B), and A s -> A S (B). 



(22) 



It is found that deviations of bp and 6^ from those determined from the D(D*) mesons affect the 
fs s I /b ratio only by small amounts. Changing bo and 64 by ~ 15% changes f Bs / fs ratio by only 
about 1%. On the other hand, the effect of g BB * w is much larger. If heavy quark flavor symmetry 
were exact, one would have gsB*-K = 9dd*-k = 0.6. However, lattice QCD simulations indicate 
that gBB*n is most likely smaller than goD*-K- For instance, two most recent Nf = 2 studies give 



g BB ^ = 0.516(5) (33) (28) (28) M\ and g BB * n = 0.44 ± 0.03^ |2fJ. Using 0.516 as the 



central value and 0.60 (0.44) as the upper(lower) bounds for gsB'-K, we find: 



1 99 +0.05 
i - zz -0.045 



(23) 



which is compatible with the result of Ref. 113 311 but with larger uncertainties 
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FIG. 4. (Color online) Same as Fig.[3] but on the left panel only the lattice extrapolations with m q / m s < 0.2 
are fitted while on the right panel only those with m q /m s > 0.2 are fitted. 

which agrees very well with the most precise result from the HPQCD collaboration: fs a / fs = 
1.226(26) 113711 . The uncertainty of ~ 0.05 does not take into account all sources of uncertainties^, 
but nevertheless it represents a reasonable estimate by covering a range of Qbb*-k values suggested 
by the two recent 1QCD calculations and possible 1/ttiq corrections to b D and 64. 



IV. SUMMARY AND CONCLUSIONS 



We have derived a co variant formulation of ^PT in order to study the D and D s decay constants. 
To simplify the analysis, we have taken extrapolated lattice data (the HPQCD continuum limits) 
as a benchmark and focused on the light quark mass evolution of f D and f Ds , in particular on 
the SU(3) breaking pattern. We find that covariant %PT describes the HPQCD extrapolations 
considerably better than HM^PT at a given order, although both approaches show improvement 
when going from NLO to NNLO. Our studies show once more that if the lattice simulations are 
performed with relatively large light quark masses (e.g., > 300 MeV), a covariant formulation 
of %PT is a better choice for chiral extrapolations, particularly at low chiral orders. 

Lattice QCD calculations have made remarkable progress in recent years. For "gold plated" 
physical quantities such as the decay constants studied in this work, the overall uncertainties have 
been reduced to a few percent. Chiral perturbation theory plays an important role in understand- 
ing some of the systematic errors, such as those from finite volume and extrapolations of the light 
quark masses to their physical values. On the other hand these precise lattice data are also valuable 
to fix the relevant LECs appearing in xPT, which can then be used to predict physical observables 
8 For instance, the small uncertainties propagated from the 1QCD results of fjj (JdJ- 
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involving the same LECs. In the present work, we have used the HPQCD D (D s ) data in combi- 
nation with two lattice determinations of qbb*-k in order to predict the ratio JbJ Jb = l^l";^. 
This ratio turns out to be more sensitive to the value of qbb*-k than to possible 1/mg corrections 
to the two relevant LECs, b D and b A . 
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VI. APPENDIX 



The functions A and B appearing in the calculation of the D and D s meson decay constants 
in the text are defined as, respectively: 



A (m 2 



-m 2 log 



m? 



16tt 2 
1 



16tt 2 
1 1 



m i lo g(^)- m i 



16tt 2 

if p\ = and mi = m 2 

if p 2 — and m x ^ m 2 



(24) 



log 



16tt 2 2pJ 

+2^2m\ {pj + m\) — (m% — pj) 2 — m\ 



U + (ml - ml - pi) log 



x < tan 1 



TO 2 -m 2 -p 2 



^ 2m 2 ( P 2 + m 2 )-(m 2 -p 2 ) 2 - m 4 



— tan 1 



^ 2m 2 (p 2 + m 2 )-(m 2 -p 2 ) 2 -mt 

(25) 



In the present work, the loop results are regularized using the modified minimal subtraction scheme 
and, unless otherwise specified, the regularization scale ji is set at 1 GeV. 
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